Diffusive magnetotransport in a 2D Rashba system 
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We present calculations of the conductivity tensor a of a 2_D-system with the Rashba spin-orbit 
interaction (SOI) in an orthogonal magnetic field, with allowance for electron elastic scattering by a 
Gaussian (5-correlated random potential in the self-consistent Born approximation. The calculations 
are performed proceeding from the Kubo formula using a new exact representation of the one- 
particle Green function of the 2_D-system with SOI in an arbitrary magnetic field. We have obtained 
the analytical expressions for the density of states and a which have a simple interpretation in terms 
of the two-subband model and hold good in a wide range from the classical magnetic fields {uJc'^ <C 1) 
up to the quantizing ones (WcT ^ !)• The numerical analysis of the Shubnikov — de Haas oscillations 
of the kinetic coeffitients and of their behavior in the classical fields region is performed. 

PACS numbers: 71.10.Ej, 72.15.Gd, 73.20.At, 73.21.-b 



I. INTRODUCTION 

The growing interest in studying the spin-orbit inter- 
action (SOI) in semiconductor 2r'-structures is mostly 
due to its potential application to the spin-based elec- 
tronic devicesi. There are two main types of SOI in 
the quantum well based on zinc-blende-lattice semicon- 
ductors. First, the Dresselhaus interaction^ that orig- 
inates from the bulk inversion asymmetry (BIA); sec- 
ond, the Rashba interaction'^ induced by structural inver- 
sion asymmetry (SIA) of the confined field of a quantum 
well. Both these interactions lead to the momentum- 
dependent spin splitting of the electron energy spec- 
trum and to the formation of quantum states with 
the hard linked spatial and spin degrees of freedom of 
the electrons. They are responsible for many interest- 
ing effects in the transport phenomena like beatings in 
the Shubnikov — de Haas (SdH) oscillations^ii; weak 
antilocalization^*S*2iS; current-induced non-equilibrium 
spin polarization^iii; spin Hall effec1?i2ii^, and so on. 

At present there are some sufficiently well developed 
theories of the kinetic and spin phenomena in 2D- 
systems with SOI in zero or classical weak {uj^t ^ 1) 
orthogonal magnetic fields. Here uj^ — \e\B/mc is the cy- 
clotron frequency, and r is the electron scattering time. 
As for theoretical studies of the considered systems in 
strong, and especially in quantizing {uj^t > 1) magnetic 
fields, there is still no satisfactory analytical descrip- 
tion of the kinetic phenomena even in the usual diffu- 
sive regime (without quantum corrections). The com- 
plex form of the eigenspinors and energy spectrum of an 
electron in the presence of SOI and a strong magnetic 
field^i is the main cause of such a situation. Direct em- 
ployment of this basis forces one to proceed almost right 
from the start to the numerical analysis of very cumber- 
some expressionsi^iiSiiSiii. 

The strong magnetic field is however one of the most 
efficient tools for investigation of SOI* and manipulation 
of the spin degrees of freedom in the semiconductor 2D- 
structures. Thus, a rather simple theoretical description 



of the kinetic phenomena in the 2£)-systems with SOI 
in a strong orthogonal magnetic field becomes a neces- 
sity. In the present work we consider the problem of 
calculation of the longitudinal and Hall resistances of a 
disordered Rashba system in the self-consistent Born ap- 
proximation (SCBA). 

We have found the exact relation between the one- 
particle Green function (GF) of the Rashba 2r'-electron 
in an arbitrary orthogonal magnetic field and the well 
known GF of an "ideal" electron, that is an electon with 
the ideal value of the Zeeman coupling (^q — 2) and 
without SOI. This allows one to obtain the analytical 
expressions for the density of states (DOS) and the con- 
ductivity tensor a^j in the SCBA. These expressions hold 
good in a wide range, from the classically weak magnetic 
fields {iv^T ^ 1) up to the quantizing ones {uj^t > 1). 
They have a simple physical interpretation in the frame- 
work of the two-subband conductor model. On the basis 
of these results, we perform the numerical analysis of the 
beatings of the SdH oscillations of the considered kinetic 
coefficients, as well as of their behavior in the classical 
magnetic fields region. 



II. MODEL 

Let us consider a two-dimensional (||OXF) degener- 
ate gas of electrons with effective mass m, and effective 
Zeeman coupling g that moves in a Gaussian (5-correlated 
random field U{r) in the presence of an external orthog- 
onal (B||OZ) magnetic field B = V x A. We assume the 
Rashba interaction to be the dominant mechanism of the 
energy spin splitting in the absence of a magnetic field. 
This situation occurs, for example, in the narrow-gap 
semiconductor heterostructures, such as InAs/GaSl^, 
InGa As/In Al As The one-particle Hamiltonian of 
the considered system has the form 



H + U 



1 



2m 



+ a{(T X tt) • n + -guj^a^ + U (r) (1) 



2 



{H = 1). Here tt = p — eA/c = mv is the operator of 
the kinematic electron momentum; cr = (cr^, a^) is 
the vector formed by the Pauh spin matrices; a is the 
Rashba spin-orbit couphng; g is the effective Zeeman 
couphng (f;~factor). 

In the gauge A = {0,Bx,0), the components of the 
eigenspinors of the Hamihonian Q of a free {U{r) = 0) 
Rashba electron are expressed through the Landau wave 
functions V'„ x i^) depending on the Landau level number 
n — 0,1,2,... and the X-coordinate of the cyclotron 
orbit centre X = —ky/muj^ 



that is diagonal in the basis (|2a|l and approaches the 
helicity operator (cr x p) • n/|p| as i? ^ 0. Here n is the 
unit normal vector to the considered 2D-system; 



^ 1 



(4) 



1 



C„„V'«-i,x(r) 



1 + L V'„,x(r) 



, a — {s, n, X) . 

(2a) 

The corresponding energy levels have the following form 



f — 



— o^ciJ, n = 0, s = +1 . 

uj^ [n + Sy/S"^ + ] , n > 0, s = ±1 . 



(2b) 

Here C^ .^ = ^^/2n/ [s^Jd"^ + 27^71 — 5\ is a normalizing 
coefficient; 5 = {g — 2)/A is the relative deviation of the 
effective Zeeman coupling from its ideal value ffg — ^ (for 
definiteness, it is assumed that (5 < in these equations, 
but all the following results are valid for any sign of (5); 
and, finally, 7 = a^Jm/uj^ is the dimensionless Rashba 
spin-orbit coupling. 

The quantum number s = ±1 describes the helicity 
of the Rashba electron eigenstate as in the absence of a 
magnetic field^^. Indeed, it can be verified immediately 
that s = ±1 is the eigenvalue of the operator 



acr X TV 



(3) 



is the Hamiltonian of the "ideal" {g^ = 2) electron in a 
magnetic field, which commutes with cr • n, (cr x tt) • n, 
and with H (jlj. 

In spite of this analogy with the B = case, we cannot 
say that the Rashba electron has in the states (|2a|) the 
spin projection ±1/2 onto the direction an x n + Lu^6n, 
because the components of the kinematic momentum op- 
erator TT are not commuting motion integrals. Neverthe- 
less, this interpretation makes sense in the quasiclassical 
limit, when one can speak about the electron path in a 
magnetic field. Namely, the quantum number s = ±1 de- 
termines the value of the spin projection on the instant 
direction of aTT x n -f oj^Sn that changes along the quasi- 
classical electron path. Thus, the spin configurations of 
the Rashba electron states form vortices in the XY-p\a,ne 
with center at the origin. 

The conductivity tensor a of the considered system 
has just one independent circularly polarized component 
(T = cr^^ -|- iay^. In the one-electron approximation, it 
has the formic 



-a^^^^Try, 



[2<i>f ^ - 



AA] I 

EE\ \e=E^ 



{dE-dE-)[<^rE'-^%%' 



dE 



(5) 



E'=E 



^EE' 



{G^ {E)V_G^ [E')) is the current ver- 



Here, 

tex operator; = zL iVy — Vj. zL 2iaaj. are cir- 
cularly polarized components of the full velocity oper- 
ator [the corresponding components cr are defined as 
CT_j_ = (cr^ ± iay)/2], where the last term occurs due to 

SOI CEJ. G^^^\E) = l/{E~n-U±iO) is the resolvent 
(retarded {R) or advanced {A)) of the Hamiltonian 
and angular brackets (. . .) denote the averaging over the 
random field U configurations. 



III. ONE-ELECTRON GREEN FUNCTION 

By definition, the one-particle GF is the averaged re- 
solvent of the Hamiltonian (d {G^'^^^E)) = {1/{E - 
H ~ U ± iO)). It is connected with the electron self- 
energy operator Y.^'^^') [E) by the relation {X ~ R, A) 



{G'^iE)) 



(G^iE)) 



(Gf^iE)) (Gf^iE)) 



E-n 



T,^{E) 
(6) 

The direct employment of the eigenspinors Q for cal- 
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culation of JB)), or kinetic and thermodynamic proper- 
ties of the Rashba system in a strong magnetic field 
leads to very complicated expressions. One is forced 
almost from the first steps either to turn to numeri- 
cal calculationai^iiSiii, or to make simplifying approxi- 
mations like the momentum-independent spin-splitting 
energyi^. This makes more difficult the interpretation 
of the results obtained in such a way, as well as the un- 
derstanding of the whole physical picture. But it turns 
out that the GF of the free (C/ = 0) Rashba system is 
expressed exactly through the GF of the "ideal" electron 
in a magnetic field. This opens up new possibilities for 
analytical studies of the considered system. Indeed, it is 
easy to check that the Hamiltonian of the free Rashba 
systems can be presented in the following form 



(7) 



Here v is the helicity operator defined in Eq. |j3Jl, Tig is 
the Hamiltonian of the "ideal" electron 

The substitution of the Hamiltonian lO into the resol- 
vent G[E) = {E — H)~^ gives, after some simple algebra, 
the following result (here and below, we drop superscripts 
R{A), if this does not lead to misunderstundings. Some- 
times, for brevity of notations, we shall not write explic- 
itly the energy arguments of the resolvents or GF's.) 



G{E) 



E -Hq+ [a(7r X n) -f tt^c^"^] ' 



[E -f ma^ 



■Ho)' 



(8) 



the free Rashba system 
GiE) = 

1 ^ f2^ -f 4s[ma^ — w^Jcr^ — a(7r X n) • <t] 



2n 



s=±l/2 



E + ma^ + sVlg - 



E 

=±1/2 



aiiT X n) • (T 



G{E + ma' + sVlg) . 

(10) 



We use here the same notation (G) for the GF of the 
Rashba electron and for the GF of the "ideal" electron. 
However, this does not lead to confusion since the latter 
depends always on the energy arguments like E + ma' + 
sftg etc. 

It is important that the same representation can be ob- 
tained for the averaged resolvent of the Rashba system in 
the SCBA. We restrict ourselves here to an approxima- 
tion in which the electron self-energy operator is diagonal 
in the spin space. Then, the SCBA equation for Y.^ {E) 
has the following form 



t{E) = W{SpG{E)) = 



Y^^{E) 








(11) 



where Sp denotes the trace only over the spatial de- 
grees of freedom; W is the mean-square fluctuation of 
the Gaussian random field ((C/(r)[/(r')) = WS{r — r')). 
Therefore, it suffices to make everywhere in Eq. (|10|) the 
following substitutions 



E ^ E-Y,iE) 



g^c 



g^c 



-4S„(£;) (12) 



where 

= 2y^2ma^E + m^a'^ -\- u^S^ = ^Vt^ + AuolS"^ . (9) 

The quantity is equal to the magnetic field- 

dependent frequency of the spin precession of the elec- 
tron with energy E that is responsible for the Dyakonov 
— Perel spin relaxation mechanism.^''; Vl is the same fre- 
quency in the absence of a magnetic field. It should 
be noted that the same representation of the one- 
electron GF can be also obtained for a system with the 
momentum-linear Dresselhaus SOI. For example, in the 
case of a [001]-grown quantum well based on the AjjjBy 
semiconductors, it suffices to replace tt ^ tt = (7ry,7r^) 
in the definition of the helicity operator ©, change the 
sign before the Zeeman term {g^ = — 2!) in the Hamilto- 
nian of the " ideal" electron Q), and, finally, to redefine 
the parameter 5 ^ 5jj = {g + 2)/4. 

The denominator of the right-hand side of Eq. ((SJ de- 
pends on the "ideal" electron Hamiltonian alone. Ex- 
panding this expression into the partial fractions, we ob- 
tain the desired representation of the one-electron GF of 



to obtain the desired representations for the averaged 
GF's in the SCBA. Here E^(^) [E) = [S^ [E) ± {E)] /2 
are the even and odd parts of the electron self-energy. 
The first (E^ = ± i/2T^) describes the perturbation 
(shift Ag and broadening I/t^) of the one-electron energy 
levels by a random field. The real part of E^ = /S.^±i/2T^ 
defines the renormalization of the Zeeman coupling H12|l . 
while its imaginary part cx 1 /t^ makes a contribution to 
the overall broadening of the one-electron energy levels. 
As a result, we obtain a expression like Eq. (|10|l for the 
averaged GF, where 

QRiA) + jna' + sns) = ^ - 

E + ma^ + sVLr, -Ho±^ 

(13) 

is the averaged retarded (advanced) GF of the "ideal" 
electron, and 
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are the disorder-modified frequency of the spin preces- 
sion and the inverse Hfe time of an electron in the s-th 
spin-split subband. As usual, we do not take explicitly 
into consideration in the one-electron energy levels 
shift Ag that is absorbed by the normalization condition, 
but we mean here that the odd shift is included in 
the definition of the effective (7-factor in accordance with 
(|12|l . The explicit allowance for the Zeeman coupling 
renormalization is particularly important in the SdH os- 
cillations regime. 



IV. DENSITY OF STATES AND SELF-ENERGY 

We first consider the calculation of the DOS 'n{E) = 
Im(Tr G^{E)) /it using the above-obtained expression for 
the one-particle GF Here, the symbol Tr denotes 

the trace over the spatial and spin degrees of freedom. 
For the sake of simplicity, we shall deal with the case of 
large filling numbers [E ^ ujj. Calculating the trace of 
resolvent (|10|l over the spatial and spin degrees of free- 
dom, we obtain the following expression for the DOS 



n{E)^ —n^°^[E + ma^ + sinB±^c)] 

s=±l/2 

= Y: "^nfKE). (15) 

s=±l/2 

Here, we take into account that the DOS of a spinless 
electron in an orthogonal magnetic field rt(°^ [E) satisfies 
n(o)(£;) = n^°\E ± cjj at large filling factors {E > wj. 
The sign before ujc is chosen in such a way as to ensure 
the right-hand limit s(r2g ztw^) — > ±s(7u;j,/2, as the spin- 
orbit coupling approaches zero. The effective mass rrig in 
the s-th subband is defined as 



m 1 



4TOa^ 



mdEiE + sflB)- (16) 



In the considered case this expression coincides with the 
usual definition of the transport and cyclotron effective 
masses in the isotropic nonparabolic bandSi. The symbol 
denotes the derivative with respect to energy E. 
In full accordance with the two-subband model, the 
DOS in Eq. H15|) is presented as a sum of partial con- 
tributions. Using this expression for the DOS, we can 
obtain the analytical form of the equation for the elec- 
tron concentration n — /^^ n{E)dE that is the normal- 
ization condition for the Fermi level determination. For 
example, at _B = we have 



m 

TT 



Ep + ma j — — Eq 



(17) 



Thus, the energy i?Q ~ Ep 



corresponds to the 



Fermi level in the absence of SOI. Notice that the partial 
electron concentrations = rn{EQ + sO^)/27r depend 
nonlinearly on the Fermi energy, in contrast to n (|17|l . 
Of course, the difference between Eq and Ep is small for 



weak SOI (ma^ <C Ep), but it should be taken into ac- 
count when analyzing the SdH oscillations that are very 
sensitive to the electron spectrum character. 

The representation H15|) allows one to obtain a simple 
analytical experssion for the DOS that holds good up to 
the quantizing fields region {lli^t > 1). Indeed, the DOS 
of a spinless electron in the large filling factors region 
{E ^ oj^) has the form 



sinh ■ 



n(")(^) 



2tt ,71- 
cosh 



cos 27r- 



E 



(18) 



Inserting Eq. (|18() into Eq. H15() . we obtain for the oscil- 
lating part of the DOS the following expression 



* , , 2m 
An(Epj = exp 

TT 



TT 
UJ^T 



2ma 



Er) r2 f, 

cos27r— ^ cosTT — —- 
■ sm ZTT — - sm TT —2 



(19) 



that is valid in the magnetic fields region under consider- 
ation. It follows from ifT^I) that energy Eq (see Eq. ^7^) 
defines the main period of the SdH oscillations, and 
(see Eq. Q) defines their beating period that depends 
on the magnetic field (see Fig. 1(a)). The second term in 
Eq. H19|l appears due to the difference between the effec- 
tive masses (|16|l . In the case of weak SOI (fl <C E), 
the oscillations of the DOS are determined completely by 
the first term in Eq. H19|l . Then, the location of the fc-th 
node of beatings is determined by the condition 



2n 



v/(2fc + 1)2 -(.9-2)2 



(20) 



This limit was considered in Ref'^^. Unlike the results 
of that work, the above-obtained equations still stand in 
the case of strong SOI. In addition, we have taken into ac- 
count the Zeeman splitting of the electron spectrum that 
allows to describe more correctly the oscillation pattern. 
For example, the Eq. (|20|l allows to determine both the 
spin-orbit a and Zeeman g couplings by measured loca- 
tions of two different nodes (see upper curve in Fig. 1(a)). 
On the other hand, the spin precession frequency fig ap- 
proaches \6\uj^ as the magnetic field B increases. There- 
fore, in this case a gradual transition from the beatings 
of the SdH oscillations to the familiar Zeeman splitting 
of the oscillating peaks should be observed. The begin- 
ning of this transition can be seen on the lower curve in 
Fig. 1(a). 

Another important characteristic of the one-electron 
states of the 2_D-Rashba system is the difference of the 
partial DOS's with opposite spin projections onto the 
OZ-axis 

6niE)=n^,{E)-n^iE)^^ J2 ^^1"^^?) (21) 



s=±l/2 
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(as 



FIG. 1: Plots of the SdH oscillations of the total DOS (up- 
per panel) and of the difference of the partial DOS's (bot- 
tom panel) of the 2D Rashba system at fixed g = 2.8, and 
kpl — 35.0, and different fir = 3.0; 1.5; 0.75 (up to down). 
The arrows point the nodes location that are calculated with 
Eq. CT . 




FIG. 2: Plots of the SdH oscillations of the inverse life time 
of one-electron states in the s-th spin-splitted subband at 
different values of Zeeman factor g = 1.8; 1.0; 0.2 (up to 
down), and fixed values of kpl = 35.0, and fir = 1.5. 



This quantity is proportional to the derivative of the 
transverse spin magnetization with respect to energy E 
and, therefore, it enters in the dcfinitition of the effective 
concentrations of current carriers in the dissipative part 
of the 2Z3-Rashba system conductivity in an orthogonal 



magnetic field (see the next section). 

Evidently, the difference of the partial DOS's (|21|l van- 
ishes in the region of classical magnetic fields {oj^t <C 1), 
but it plays an important role in the SdH oscillations 
regime. In the case of large filling fsctors, the oscillating 
behavior of this quantity is described by the following 
expression 



2m2uj.S / TT \ Er. n„ 
— — exp sin 2tt — sin vr — — . 

(22) 

Unlike the total DOS 119|l . this expression contains just 
one oscillating term, because 6n(E) does not depend on 
the effective masses llt)|) . Indeed, the difference of 
the partial DOS's Sn{E) is non zero, which is entirely 
due to the spin degrees of freedom of the electrons. The 
typical SdH oscillation patterns of Sn{E) are depicted in 
Fig. 1(b). 

Now, let us turn to the discussion of the electron life 
time in the s-th spin-split subband which is defined, 
according to Eq. 1)14(1 , by the imaginary parts of the even 
and odd self-energies ^e(^o}- other words, the total life 
time of the one-electron states is determined by the 
sum of the weighted relaxation rates of the orbital and 
spin degrees of freedom. The first term in this expression 
is proportional to the above-considered total DOS, hence 
its magnetic-field dependence coincides up to the scale 
factor with the patterns shown in Fig. 1(a). Of particular 
interest is the last term in Eq. (|14|l stemming from the 
Zeeman coupling renormalization. It is proportional to 
the difference of the partial DOS's and, therefore, 
plays an important role in the SdH oscillation regime, as 
shown in Fig. 2. Notice that the beatings of the SdH 
oscillations get supressed with the increase of the rela- 
tive magnitude of the second term in Eq. H14|) . Indeed, 
Eq. H2l)|) determines the location of the beating loops of 
the oscillation instead of the nodes. Thus, the broaden- 
ing of the Zeeman levels leads to observable supression 
of the beatings of the SdH oscillations, as in the case 
of the competition between the Rashba and Dresselhaus 



V. CONDUCTIVITY 

The general expression for the conductivity (jSJ con- 
sists of two different terms. The first of them describes 
the contribution of the electrons at the Fermi level, the 
second one contains the contribitions of all filled states 
below the Fermi level. We begin the calculation of the 
conductivity with the last term of |(5j cr^^. First of all, it 
is pure imaginary and, therefore, makes a contribution in 
the Hall conductivity alone. Stfeda#i was first to show 
that, in the absence of SOI, this part of the conductivity 
is equal to 



a" =i\e\c 



dn 



(23) 
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where n is the electron concentration. It should be 
pointed out that Eq. (|23|) is exact, and with the thermo- 
dynamic Maxwell relation <t^^ can be expressed through 
{dM/dE)g, where M is the orbital magnetization of the 
electron gas. Detailed discussion of a^^ and its physical 
interpretation can be found in surve}*^^. 

This result is extended immediately to the electron sys- 
tems with SOI. Following Sti^eda's argument, it can be 
shown that the part a^' of the of 2D Rashba system 
conductivity is expressed as 



a''=i\e\c 



-i\e\c 



[ dE 



dn 

dn \ _ g\e\ 
dB ) „ imc 



(s = ±1/2) are direct analogues of the familiar param- 
eter = En^'^''{E) that stands for the current carrier 
concentration in the dissipative part of the conductivity 
tensor of spinlcss 2L'~electrons in a magnetic field in the 
SCBA-'^^. In the classical magnetic fields region, they are 
equal to the partial electron concentrations in the spin- 
split bands — ttl^Eq -\- sflg)/2'!T. On the other hand, 
(EH) approaches {E - sguj^/2)n'^°\E - sguj^/2) in the 
limit a — > 0, in which only the Zeeman energy splitting 
remains. The first two equalities in Eq. H24|l are also ex- 
act. The last one is obtained in the SCBA at large filling 
factors {E ^ u>^) using the expression ((TB|l for the DOS 
of the 2D Rashba system in a magnetic field (for details 
see Appendix B) . It should be emphasized that a similar 
SCBA expression is valid also for spinless electronaSi. 



.\e-\ncv 



(24) 



where is the spin magnetization of the electron gas. 
The quantities 



2^ 



n(o)[i?o + .s(f7s-c.J] 
(25) 



Now, we turn to the first term in the conductivity |(SJ). 
It is quite easy to show, by identical transformations, 
that 



47r 



ReTrV^^ii ^ _^ReTr(G^) = ^^ReTr \<^>f<ffG^iE, + sns) 



27rm 



s^s) , (26) 



where the averaged GF's of the "ideal" electron are de- 
fined in Eqs. H13() and H14|l . In obtaining the last equality 
in H2t)|) . we used the immediately verified identities 



(27) 

The main contribution to the dissipative part of the con- 
ductivity is proportional to the current vertex in 
Eq. ©. If we accept the SCBA I^IJ for the electron self- 
energy E, we must evaluate the conductivity © in the 
ladder approximation in order to satisfy the particle con- 
servation law. But as show the calculations, the relative 
magnitude of the ladder correction to the conductivity is 
second-order in the small parameter l/{kpl), and it can 
be neglected in comparison with the "bare" part of the 
conductivity (Acr'^^'^/cr = 0.01 0.001 for typical values 
of kpl = 10^30). 

Thus, it suffices to calculate the "bare" part of the con- 
ductivity that is obtained by replacement TvVj^'^^j^ 
Try+(G^)y_(G^) in the first term in the right-hand 
side of Eq. (|SJ . We drop the details of calculations that 
can be found in Appendix B and proceed to the results. 
The overall contribution of the three above-considered 



parts has the usual Drude — Boltzmann form 



B 



i^i^B 



(28) 



(it is meant here and below that B B/c). Here, the 
first two terms represent the sum of the partial conduc- 
tivities of the electrons of two subbands with different 
mobilities /^^ ~ \e\Tg/m, and effective concentrations 



TV, = TV, + 



rs^B (Stt) 



;{E - sns)^ilB^ 



(29) 



This expression differs from (see Eq. (|25|) ') by the 
second term that originates from the principal values of 
the one-electron GF's. 

The last term in Eq. H28II represents the small correc- 
tion to the conductivity due to the electron intersubband 
transitions. In the leading approximation in powers of 
the smallness parameters lu^/E and fl^/E, it has the 
form 



Act 



lei 2ma^ 



B 



n2 



-U! 



l(0) 



(30) 
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where n(°) = J2s "-s"^ (E) is the DOS at the Fermi level 
averaged over the electron subbands. The relative con- 
tribution of this correction to the full conductivity l|28|) is 
of the order of magnitude (uj^/E)^ and can be neglected 
in the large filling factors {E ^ u;^) region. 

We emphasize that the contribution to the conductiv- 
ity from intersubband transitions vanishes as the mag- 
netic field approaches zero. This would be expected, be- 
cause the conductivity tensor in the absence of a mag- 
netic field is diagonal in the original spin space {a^^ = 
= 0) by virtue of the momentum parity of the GF's, 
and the full conductivity is equal to cr = a-^.^ + ctjj^^^- In 
fact, the case in point concerns the time inversion sym- 
metry. Using a unitary transformation, it can be turned 
into a matrix s-representation in which the one-electron 
GF's are diagonal and, therefore, a — Cji+Cj^ — cr_^_ + a_ 
due to the trace conservation. Applying an external mag- 
netic field breaks the above-mentioned symmetry that 
is responsible for the appearance of the intersubband 
transition-induced conductivity Act. 



VI. RESULTS AND DISCUSSION 

First of all, let us summarize briefly the main results 
obtained in this work. We have shown that the eigen- 
states of the 2D Rashba electron in an orthogonal mag- 
netic field are characterized by a special motion integral 
(0) that generalizes the notion of helicity^. Using this 
fact, we have found the relation 1)101) between the GF's of 
the 2D Rashba electron and the "ideal" one that holds 
good for arbitrary orthogonal magnetic fields as well as 
for the strong spin-orbit coupling. With the help of this 
relation, we have obtained, in contrast to Refs.i^ji^tii, 
the analytical SCBA expressions for the DOS H15|) and 
magnetoconductivity (|28|l of the 2D Rashba system that 
are valid in a wide range from the classical magnetic fields 
up to the quantizing ones {lo^t > 1). They permit a sim- 
ple interpretation in the framework of the model of two 
types of current carriers with different concentrations and 
mobilities. The spin-orbit as well as the Zeeman splitting 
of the electron energy are properly allowed for in these 
expressions, unlike the results of Ref.— . We have shown 
that the competition of the relaxation rates of the or- 
bital and spin degrees of freedom in the total inverse life 
time 1 /Tg of the one-electron states in the s-th subband 
leads to the supression of beatings of the SdH oscillations 
as does the competition of the Rashba and Dresselhaus 
SOTsi^. Finally, we have shown that the breaking of the 
time inversion symmetry in a magnetic field leads to the 
appearance of the intersubband term in the 2D Rashba 
system conductivity if^ . 

We start the discussion of the results with the con- 
ductivity in zero magnetic field. In this case, it follows 
immediately from 1)28(1 that 




|e|(n, /i, 



1 



/ mi 



(31) 



FIG. 3: The smooth magnetic-field dependence of the re- 
sistance calculated with Egs. 1281 — 1301 at fixed g = 3.5, and 
kpl = 32, and different Q,t = 1.5; 1.0; 0.5 (up to down). 



Thus, the Rashba spin-orbit interaction leads to a de- 
crease in conductivity, and not to its increase, as it was 
claimed in Ref.^^. Let us note that the authors of that 
work ignored the difference in mobility between the elec- 
trons of different subbands. As a result, thay obtained 
a correction to the conductivity of opposite sign com- 
pared to ()31|l . which is actually unobservable, because it 
is absorbed by normalization condition H17() . 

Let us proceed now to the discussion of the magne- 
totransport in the 2D Rashba system. It is well known 
that in the two-subband conductors the classical positive 
magnetoresistance and the magnetic-field-dependent 
Hall coefhcient are observed (see, for example, ^^). The 
considered system differs from a classical two-subband 
conductor in two points. Firstly, the mobilities and ef- 
fective concentrations of current carriers l|25|l depend on 
the magnetic field. Secondly, the full conductivity of the 
2D Rashba system in an orthogonal magnetic field is not 
an additive sum of the intrasubband contributions, but 
it contains a nonadditive intersubband term (|28|l . How- 
ever, all these factors lead to very slight magnetic field 
dependences of kinetic coefficients due to small differ- 
ences between concentrations and mobilities of current 
carriers. For example, the relative magnitude of the clas- 
sical magnetoresistive effect comes to only 1-^2% for 
typical values of parameters (see Fig. 3). 

In discussing the SdH oscillations, we restrict ourselves 
to the consideration of the large filling factors {E 3> lo^ 
region, where the SCBA is applicable to the description 
of the one-electron states and kinetic phenomena. As 
usual, we extract in the linear approximation the oscil- 
lating parts of the conductivity that enter through par- 
tial DOS's into the effective concentrations (|25|l and 
mobilities /i^. Neglecting the small differences between 
concentrations and mobilities of current carriers in the 
smooth parts of conductivity, we obtain the expressions 
for the oscillating parts of the longitudinal p resistance 
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FIG. 4: Plots of the SdH oscillations of the longitudinal 
magnetoresistance (upper panel) and Hall coefficient (bot- 
tom panel) of the 2D Rashba system calculated with Eqs. 1281 
— I3(JI at fixed g = 3.5, and kpl = 32, and different fir — 
1.5; 1.0; 0.5 (up to down). 



and Hall coefficient R 



Ap(i3) 
P(0) 
AR 



H 



:4exp 
2 



H 



TT 
UJ^.T 



COS 27r 



-^COSTT— ^ 



(32a) 



R%{0) n^B^ 



cxp 



TT 
UJ^.T 



COS 2TT 



En ^ D 

-^COSTT— ^ 



(32b) 



Here, p(0) = l/a'(0) and R^ ~ —l/\e\nc are is the resis- 
tance (see Eq. (|31|) ) and Hall coefficient in zero magnetic 
field respectively. Results of numerical calculations of 
the SdH oscillations are performed using total expessions 
— (1201) are showed in Fig. 4. 

Up to definition of the beatings period, the Eqs. H32fl 
agree with expressions for the longitudinal and Hall con- 
ductivities obtained in Ref.^^. As pointed out above, the 
spin precession frequency Qg and, therefore, beatings 
period (see, for example, O) depend on the magnetic 
field due to allowing for Zeeman coupling. As result, the 
measurement of two different nodes location permits to 
determine both the spin-orbit a and Zeeman g couplings 
using the Eq. ((201 . 
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APPENDIX A: SOME USEFUL IDENTITIES 

We obtain here several identities for the off-diagonal 
matrix elements (G^j^q^-j) of the one-particle GF ^ that 
are necessary for calculation of the kinetic coefficients 
(see, for example, Eq. (|B5I) ). 

The matrix of the one-electron GF satisfies the Dyson 
equation 

^2 1 ] . „ . 

(G) - an - (o- X 7r)(G) =/, 



2m 



(Al) 

where / is the unit 2 x 2-matrix. It is assumed that the 
odd part of the electron self-energy is included into 
the effective Zeeman coupling g (see (I12II ). The tf matrix 
element of Eq. (|Aip has the form 



E 



TV 

2m 
2m 



(G^)+za7r+(G^^> = 1.(A2) 



Analogous relations can be obtained from the conjugated 
Dyson equation (|A1(I . From their comparison with (|A2|I 
it follows that 

^-(Gxt) = -{Gn)^+ : ^+(Gtx) - . (A3) 

Now, we replace in (jA2|l the matrix elements {Gy^f^i^yi 
with their expressions through the one-particle GF of 
a "ideal" electron (see Eq. (I10|l ). As a result, we have 
the following useful relations between (G|m|p and the 
one-particle GF of the "ideal" electron 



s=±l/2 



'-'TT(ii)(^o+*^i 



, (A4) 

By combining these identities with the corresponding di- 
agonal matrix elements of Eq. (|10|l . we obtain 



'^T(^iT(Ti))T2«ma(G-f^(^)) 



±4dma ■ 



s=±l/2 



{E, + sQs) (A5) 



APPENDIX B: CALCULATION OF 
CONDUCTIVITY 

We first consider the derivation of the last equality in 
Eq. (|24|l . The immediate differentiation of the electron 
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concentration with respect to the magnetic field gives the 
following result 



n 
B 



(Bl) 



Now we should calculate the last term in this expression. 
Using the representations ifTTijl . ifT^ for the averaged GF, 
we write down 



s=±l/2 



1 



+4s- 



ma 



a{(T X tt) • n — cJ^Sa^ 



m 

'dB 



, (B2) 



Let us multiply together the expression in square brack- 
ets and the derivative of Hamiltonian H Q , (0) , keeping 
the terms diagonal in the spin space and neglecting the 
terms linear in a^. These latter make contributions pro- 
portional to ■nP^\E) — rS'^\E±ijj^ and vanish in the mag- 
netic field region of interest to us. As a result, Eq. ljB2|) 
takes the form 



Tr(G) 



m 

'dB 



= ±1/2 



dUo 1 4s 
'dW ~ BOZ 



2 r2 
UJ^O 



(B3) 



where Hq is the Hamiltonian of the "ideal" electron |0J. 
According to the well known theorem of quantum me- 
chanics, there is the identity (9i?/9A)„„ = dE^/dX, 
where E^ is the n~th eigenvalue of the Hcrmitian op- 
erator H. This make it possible to perform the following 
substitution dn^/dB -> Hq/B in Eq. Then, on 

simple rearrangements, Eq. ljBl|) takes the form 



1 



dn \ n 
dBJ ^^B~ 2^ 

^ ^ s=±l/2 



E 



X ImTr 



Er, -\- sn. 



4s 



2 r2 



G^{Eo + sns). 

(B4) 



Eqs. 124|) . (|25|l are derived immediately from this equa- 
tion. 

Now, let us proceed to the calculation of the dissi- 
pative part of the conductivity that is proportional to 
TtV+{G^)V_{G'^) in the SCBA. Performing the trace 
over the spin degrees of freedom and taking into account 
the relations (IA3II. we write it down in the form 



+2K(Gfi) + 2^ma(G5)) {n_{Gf^) - 2^ma{G^)) 



(B5) 



Two last terms in the right-hand side of this equation are calculated using identities H10() . IjASp . We present a more 
detaled calculation of one of the two first terms in HB5() . 

For example, let us substitute, in the first term in Eq. IjBSp . the diagonal matrix elements (tf) of (|10|l for (G), and 
use the identity 



G^{Ea + sng)7r_Gi{Ea + s'^^) 



(B6) 



Then, after some simple but cumbersome algebra, we can rewrite the contribution of this term to the conductivity 
in the following form 



Trm ^ — ^ " 

s 

-im ^ ^s,TT^-s,TT 



sJl , 



1 - iuj^T^ 



■ Im Sp Gf} (£^0 + sf^s) - i Re Sp G^ {Eq + sflg)] 



2TTim 



E, + ^l{n^s-n^)-^ 

1 - i{u;^ + 2sng)T^ 



Sp [G^ {E, -sQs)- G^ {E„ + sQs) 



+ -Sp [G^(£;„ - sQs) + G^{E, + sn^)] } • (B7) 
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The last terms in curly brackets are cancelled exactly by 
the corresponding terms from Eq. (|26|) . The contribution 
of the second term from Eq. I)B5|I can be transformed in 
a similar way. 



Of course, it is necessary to perform a set of unwieldy 
some transformations to obtain Eqs. (|28|l . (|29l) . and l|30|l . 
However, further calculations have a purely technical 
character and we omit them. 
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